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\ We investigate the possibility to measure the nonforward gluon distribution function 

by means of diffractively produced 7r + 7r~ and K + K~ pairs in polarized lepton 
nucleon scattering. The resulting cross sections are small and are dominated by the 
gluonic contribution. We find relatively large spin asymmetries, both for tt + tt~ and 
for K + K~ pairs. 
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Ph. 1 Introduction 



QCD parton distribution functions are specified by matrix-elements of nonlo- 
cal quark and gluon twist- 2 operators separated by a light like distance and 
taken inside a forward nucleon matrix element. Although these forward ma- 
trix elements were generalized long time ago to nonforward parton distribution 
functions [1], where nonlocal QCD operators are sandwiched between states of 
unequal momenta, they became only recently the object of intense studies [2- 
4] . This interest was driven by the hope that the measurement of nonforward 
parton distributions (NFPDs) and subsequent extrapolation to the forward 
limit might shed more light on the spin and angular momentum composition 
of the nucleon [5]. 

Contrary to ordinary parton distribution functions, nonforward parton distri- 
bution functions are probed when the nucleon recoils elastically. Thus NFPDs 
provide a natural interpolation between ordinary parton distribution functions 
and elastic form factors. It was shown that NFPDs are in principle measurable 
in deeply virtual Compton scattering or through exclusive production of light 
vector mesons. 



TPR-99-02 



16 April 1998 



Fig. 1. The process of the quark antiquark pair production to lowest order QED. 

In this paper we shall consider the diffractive leptoproduction of quark an- 
tiquark pairs in a kinematics where the produced quark pair does not form 
a single meson but fragments into a 7r + 7r~ or K + K~ pair with compara- 
tively large pr (Fig. 1). We will consider the case of longitudinal polarized 
electron beam and proton target. This process is quite similar to diffractive 
leptoproduction of dijets, analyzed in [6,7], which was proposed to allow for a 
measurement of the polarized glue in the nucleon. We notice that at medium 
energies like in the HERMES experiment dijets cannot be observed. In this 
experiment one can only look for single hadrons representing the jets. The idea 
to access the polarized forward gluon distribution function of the nucleon via 
the measurement of a double hadronic high momentum final state has been 
proposed in [8]. 

We treat the fragmentation process of the two produced quarks into the final 
hadrons using the LUND string fragmentation model [9,10]. To be precise we 
use this model to extrapolate from dijet production, where it is known to be 
valid, to the case where the two jets are represented by e.g. n + and n~ . While 
we do not expect the LUND model to give precise results in this few particle 
process it should allow for an order of magnitude estimate. 

In our calculation we will assume that two pion (kaon) leptoproduction can 
be factorized into three parts. We have two soft parts: one describing the frag- 
mentation of the two partons, the other the production of two partons from 
the elastic scattering off the nucleon. The second soft part will be given by 
a nucleon matrix element of a light cone operator parameterized by NFPDs. 
The hard part connects the elastic scattering process and subsequent fragmen- 
tation. 

Technically, the process we consider is closely related to the exclusive pro- 
duction of a light vector meson. In such a case the proof of factorization has 
been given in [11]. Recall that the proof only holds for longitudinal photons. 
Meson production by a transversely polarized photon is plagued by endpoint 
singularities as noted in [12]. Thus a general factorization theorem does not 
exist in this case. Further the factorization theorem was proven only on the 
level of twist-2. However, it was noted that the proof in [11] can be extended 
to two-pion production if these are produced close to threshold, i.e. if the mass 
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Fig. 2. diagram for the production of a pion pair by the fragmentation of four quarks. 

of the produced final state is much smaller than the virtuality of the probing 
photon [13]. In such a kinematics it is useful to describe the pion pair in terms 
of a two-pion twist- 2 distribution amplitude [14]. The main difference to our 
calculation then is the treatment of the produced meson in terms of a light 
cone twist-2 meson wave function rather than in terms of the fragmentation 
model. In our calculation we do not restrict ourselves to the kinematic region 
M 2 (h + + h~) < Q 2 because we have to avoid too small invariant masses in 
the final quark state to justify the application of the LUND fragmentation 
model and for correspondingly large Q 2 the counting rates are too small for 
an experiment like HERMES. However, in our kinematics a hard scale is given 
by /i 2 = j4(1 + A/ 2( fe ( + +fe -) ) with pt being the transverse momentum of the 
produced quark in the center of mass system of the incoming nucleon and 
the virtual photon [15]. The hardness of this scale ensures the applicability of 
perturbation theory. The calculation therefore can be considered as a parton 
model description combined with the LUND fragmentation model. 

We would like to stress that the process we are considering is not the only 
one leading to two charged pions in the final state. An exclusively produced 
£>-Meson for instance will immediately decay into a pion-pair. However, the 
cuts in our numerical analysis make sure that we only count pions produced 
away from the ^-resonance. A pion pair can also be produced through the 
exclusive production of four quarks, each pair forming a mesonF]. A leading 
order diagram for this process is shown in Fig. 2. However, one may expect 
that such diagrams are kinematically suppressed as can be seen best in the 
center of mass frame of the incoming nucleon and the virtual photon: two of 
the four produced quarks are obtained by the decay of a gluon taken out of the 
nucleon and therefore will carry a fraction of the incoming proton momentum 
while the other quark antiquark pair carries the momentum of the virtual 
photon. Without at least one additional internal gluon line, which would lead 
to a suppression in as, it is not possible to rearrange the momenta such that 
two light mesons are formed. Also we did not consider the process of two gluon 
production and their fragmentation into a meson pair. However, it turns out 
that in the considered kinematical range contributions from gluon nonforward 
parton distributions dominate compared to quark nonforward distributions 
whereas only the latter lead to two gluon production in the leading order. 

1 We are indebted to M. Strikman for drawing our attention to this point. 
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Fig. 3. Kinematics in the boson-nucleon center of mass frame. 



The paper is organized as follows: First, we define the general kinematics, 
give the formula for the partonic cross section and derive the amplitudes for 
exclusive qq production. We then discuss the models used for the polarized 
and unpolarized NFPDs. Subsequently, we discuss the fragmentation process 
of the two final partons into the meson pair using PYTHIA/JETSET [10]. 



2 Cross Section for Diffractive qq Production 

As already mentioned we assume that meson pair production factorizes into 
hard quark pair production and a fragmentation part. We treat the frag- 
mentation process on the level of transition probabilities using the program 
JETSET. Therefore we have to calculate analytically only the cross section 
for the process 



The corresponding diagram to lowest order QED is shown in Fig. 1. For the 
definition of the particle momenta see also that figure. We work in the center 
of mass frame of the exchanged photon and the initial nucleon and neglect 
nucleon and lepton masses. Furthermore we take the momentum transfer to 
the nucleon to be proportional to the initial nucleon momentum P^—P'^ = (P^. 

— * 

We have chosen the momenta q and P to define the z- and — ^-direction 
respectively as explained in [16]. See also Fig. 3. 

The masses of the final quarks are neglected pj = p'j = 0. Recall that for 
exclusive production of only one particle the on-shell condition leads to the 
constraint that the skewedness parameter of the NFPD £ equals the Bjorken 
variable = —q 2 /2u, v = P ■ q [17]. This is clearly not the case for two 
produced particles. Therefore it is natural to introduce the variable x p = xb)/(, 
< x p < 1 which describes the invariant mass of the final qq state. Using 



e(Z) + N(P) - e(Z') + q(p f ) + q(p' f ) + N'(P') . 
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furthermore the variables 

Q 2 = -v 2 , y = ^,z f = ^,t = (P- P'f (2) 

and (p defined as the difference between the azimuthal angles of I and pf (see 
Fig. 3) we get for the differential cross section 

47 - a ™^-)V^ (3) 



dxBj dy d(dzfdtd(p Q A \47T / 



where is the usual leptonic tensor given by 



= \ E«a(07 m «a'(O«v(O7 , '«a(0 

1 A' 

= 2in v - iy - l v q» - ^Q 2 g^ - Xie^ 8(7 l e q a . (4) 

Summing over the spins of the final nucleon we define analogously the hadronic 
tensor as J2s' = T*T U with 



(2-K)%p f +p' f -q-CP)T, 

= J d'x (q(p f W f )N(P')\j,(x)\N(P)) . (5) 

j^(x) denotes the vector current that couples to the virtual photon. In second 
order perturbation theory we may write 



•2 2 

(2ii)%p f +p' f -q-CP)T, = i|- J d 4 x, J d 4 x 2 J d 4 x 3 e"^ 
x{q{p f )q{p' f )N{P')\T{j,{ Xl )j^ (6) 

with the color current j% and the gluon field A aa . We factorize quark and gluon 
fields in hard propagators and non perturbative matrix elements respectively 
using the parameterization of nonforward matrix elements of bilocal operators 
in terms of nonforward parton distributions [18]. We project onto the twist-2 
contributions by introducing a light-like vector q\ — (q + xb]P) 2 = so that 
we have for the quark matrix element 



{N{P')\i> f {xM f { Xl )\N{P)) = U(I*)i x U{P) 



x 



/ 



dX 



+/C— terms , 



(7) 
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{N(p')\i. ; (x 2 h^ 1 i. l (x 1 )\N(Pj} = mp'hsfiMP) 



x / iX 
+/C— terms , 
and for the gluon fields 



(8) 



(iV(P')|^(^)^(^i)|iV(P)) 
_ U(P f )^U(P) (_ P e qia + Paqi g " 

I 9gcr ~r 



X 



2Pgi 
dX 



2 (X-C + ie)(X-i£) 

+ 2P gi eaXv P- qi 
dX A£ f (X) 



e -iXP-Xi+i(X-C)P-X 2 _|_ e -iXP-X 2 +i(X-C)P-Xl 



X 



2 (X-C + i£)(X-ie) 
+/C— terms . 



e -iXP-xi+i(X-C)P-a;2 _ e -iXP-x 2 +i(X-C)P-xi 



(9) 



We have not written the gauge links explicitly. In the following we will neglect 
the /C-terms that are proportional to the momentum transfer to the nucleon. 
Applying Fierz transformation, restricting ourselves to unpolarized and lon- 
gitudinal polarized nucleons and taking into account the equations of motion 
we find 



N c 
i 

x 



/ 
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dX 



2P Ql v 2 z f (l-z f ) 



2 (X -( + ie) 2 (X -lef 



xl^ip^vip'f^x-cr+x^p, 

T F q^UiP'h^UiP) 1 
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2P gi v 2 z f {l-z } ) 
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(10) 



Au{Pf)H<P'f)[^Pf»- t p »] +u(pfh tl l5v(p' f )(l-2z f )A , (11) 
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T(i) = — 
" " A N c 



lC F q f ^U(P')i 1 U(P) 



1 



dX- 



2P gi v*z f {\ - z f ) 

[T{(X)+T{(X)}(2X-() 



(X-( + ie)(X - ie)(X - (x p + ie)(X - ((1 - x p ) - is) 
\u{p f )fv{p'f) [(*(* - C) + C%(1 - z f ))P, + Cx p (2z f - 

+u{p f ) 1 ,v{p' f ) V -(x{X - C) - C 2 x p (a; p - 2^(1 - z,))) J 

lCpqfMmjJiP) i 
4 iV c ^ 
l 



2Pgi 



x 

o 

X < W 



dX- 



v2z f(X - z f) 

f((x)-f{(x) 



(X-C + ie)(X - ie)(X - (x p + ie)(X - C(l - x p ) - is) 
(p f )Mp'f) [(C*(C ~X)(1+ x p (2 - Az f )) - C%(1 - z/))P M 
+ (2X(X-C) + C%)^ 
+«(P/)7^(P/)K X (C -X)- ( 2 xl)(2z f - 1)| 

- 75^ - ^75, 7m - 7^75, ^ - A^, P - -AF f ~} . 

(12) 



A typical diagram contributing to is shown in Fig. 4 on the left, an 
example for a diagram contributing to Tjf 1 in the same figure on the right. 
The color factors are denoted in the common way Cf = 4/3, Tp = 1/2, and 
the number of colors Nq = 3. qj is the electric charge of the produced quark. 
It is easy to see that by going to the limit of a massless final state, i.e. x p — 1, 
and contracting the amplitude with the polarization vector of a longitudinal 
polarized photon the results for exclusive meson production are reproduced. 
See for instance [17]. In such a kinematical limit Zf corresponds to the light 
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Fig. 4. Typical diagrams to the process involving two gluon emission and two quark 
emission of the scattered nucleon. 
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cone momentum fraction carried by the produced quark. The appearing soft 
divergences l/[zf(l — Zf)] then do not pose any problems for meson wave 
functions that vanish at the end points. In our case we have to stay off these 
divergences using cuts for Zf. Such cuts are demanded by the acceptance of 
the meson detectors because of the relation Zf « z w := Pp w /Pq. Note that 
upon contraction with a longitudinal polarization vector and the limit x p — 1 
all dependence on the polarized gluon distribution AQ cancels out. 

Contracting T*T U with the leptonic tensor L^ u leads us to the differential 
cross section -= — . — -rrr- in terms of the nonforward distributions. To the 

aXBj dy dQ dz f at dip 

polarized cross section only the interference term of polarized parton densities 
and unpolarized densities contribute. All other combinations contribute to 
unpolarized scattering. For the phenomenological analysis we have to perform 
the remaining integrations over X in the amplitudes. To do so we need some 
model for the nonforward parton distributions. 



3 Nonforward Parton Distributions 

The nonforward parton distributions used in the calculations of the previous 
section are given by an integral over the corresponding double distributions 

a 

W = jF(X-(y,y)dy (13) 


with a = min{^-, fen ^ ^ we wr ^ e a double distribution as a product of 
the corresponding usual forward distribution and a profile function 

F(x, y) = f{x)ir(x, y) (14) 
the latter has to fulfill the conditions 

l-x 

ir(x, y) = ir(x, 1 - x - y) , J n(x, y) dy = 1 . (15) 



The symmetry y <-> 1 — x — y was noted in [17]. As a model for the pro- 
file functions we take tt(x, y) = ^jSSir^ f° r the quark distributions and 
7r(x, y) = 30y j]E§j?^ f° r the gluon distributions [18]. For the (forward) parton 
distributions we use the (slightly simplified) NLO parameterizations of [19] in 
the unpolarized and [20] in the polarized case at the scale Ql = 4GeV 2 . From 
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these models we get the nonforward parton distributions using (13). In the fol- 
lowing we will neglect all effects of Q 2 evolution of the distribution functions. 
Special care is needed for the quark sea distribution for which the original 
parton distribution is more singular than 1/x-Qy In this case the integral (13) 
is not defined for X < (. The problem can be solved using the anti-symmetry 
of the remaining part of the corresponding integrand: 

i C 

jA((,X)F c (X)dX= J ^(C,C-X)^ c (C-X)dX 

o -l+c 

= 5 / A((,X) 

-i+C 

x (^(X)0(X) - F< (C - X)0(C - X)) dX 
1 

= / A((,X) 

C/2 

x (> f (X) - ^(C - X)0(C - X)) dX . (16) 
Here the divergences cancel in the difference ^(X) — ^(C ~~ X). 

In order to determine the cross section integrated over t we assume that the 
t-dependent nonforward distributions factorize into a nucleon form factor and 
the nonforward distribution at t — 

F C (X, t) = F C (X, O)F(f) = F c {X)F(t) . (17) 

This behavior is suggested (at least for small values of t) by the condition 
Jq ^((X) dX = F(t) [3] (see also [17]). As models for the form factors we take 




with a = 2, A = 0.84 GeV in the case of quark distributions and a — 3, 
A = 2.8 GeV in the case of gluon distributions. The quark form factor is 
taken from measurements of elastic electron nucleon scattering [21], the gluon 
form factor from QCD sum rule calculations [22]. Since these form factors 
decrease rapidly for increasing \t\ the integrals over t will be dominated by 
small \t\. Therefore it is justified to keep t = in the other contributions to 
the differential cross sections. In this approximation the t-integration can be 
carried out and leads to factors of the form J_ oc Fi(t)F 2 (t) dt where Fj are 
nucleon form factors for quark or gluon fields respectively. 
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Fig. 5. Probability of creating a tt + tt pair (solid circles) or a K + K (open circles) 
pair in a quark-anti quark fragmentation process with center of mass energy squared 
M 2 (h + + h~). 

4 Fragmentation into Mesons 



To describe the fragmentation process of the two outgoing quarks into a meson 
pair we use JETSET which is based on the LUND string fragmentation model 
[9,10]. In Fig. 5 we plot the probability to produce a pion or kaon pair from two 
quarks as a function of the cm-energy M 2 (h + + h~) = (pf+p'f) 2 = {(P + q) 2 = 
— x p ). For this plot we assumed that all quark-flavors are produced in 
equal ratios as soon as M 2 (h + + h~) exceeds the respective mass threshold. 

In order to evaluate the cross section for 7r + 7r~ and K + K~ pair production in 
lepton proton scattering the events generated by JETSET are weighted with 
our differential cross section for quark pair production dx ^ dy ^ d dlfi - Fig. 5 

shows that the center of mass energy squared should be less than 8 GeV 2 
to give sizable meson pair production. In order to avoid regions where we 
cannot trust the program JETSET any longer we further have to avoid too 
small values of M 2 (h + + h~). For our simulations we work in the proton rest 
frame with a lepton energy of 27.5 GeV and use an acceptance slightly smaller 
than the HERMES acceptance. We require the angles of all detected particles 
(electrons, pions and kaons) with respect to the electron beam axis to lie in 
the range 0.04 < tan 9 < 0.140. Furthermore, the detected electron should 
have an energy exceeding 1 GeV and the detected mesons an energy larger 
than 4 GeV. We also require both mesons to have a minimum transverse 
momentum in the laboratory frame of 0.5 GeV. Further we imposed the cut 
p\ > 1 GeV 2 (with px being the transverse momentum in the center of mass 
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Table 1 

The different contributions to the total cross sections. 





PTO^m cippf 7r~^~7T ( T)Y\\ 
UOD OCvjUi li li \ YJ YJ i 


prn^S SPPt" r\' J^C [ nh 1 

blUoo DbL L . ±\. J_ V I YJ YJ j 


a 99 


3.2 10 _1 


2.8 10 _1 


VAgAg 


4.2 10 -2 


4.1 10~ 2 




p. q i n — 2 
0.3 1U 


r Q 1 n — 2 
O.O 1U 




5.4 10 1 


4.0 10~ 2 


OAqAq 


3.3 1CT 3 


2.9 10~ 3 


<y q Aq 


7.4 10~ 3 


4.8 10~ 3 




-1.2 10" 1 


-8.8 10~ 2 


OAqAg 


-4.6 1(T 3 


-7.7 10~ 4 


VqAg+gAq 


-4.3 10~ 2 


-2.9 10~ 2 


total spin average 


2.9 KT 1 


2.8 lO" 1 


total spin asymmetry 


2.8 10~ 2 


2.9 10~ 2 



frame of the photon and the incoming nucleon) in order to keep the relevant 
scale /i 2 = + M 2 (h++h-) ) sufficiently large. For the remaining variables we 
used the cuts 0.02 < rr Bj < 0.05, 0.29 < y < 0.96, 0.1 < C < 0.3, and 0.04 < 
Zf < 0.96 where the lower cut in ( was dictated by numerical reasons. The 
cuts in Zf are part of the z — s jet scheme which we applied in our simulations. 
It is also used in the standard Monte Carlo leptoproduction programs [23] . 



5 Results and Discussion 

In Table 1 we have given the various contributions to the total cross section 
for meson pair production. Contributions from the polarized and unpolarized 
gluonic NFPD are denoted by the indices g and Ag, contributions from the 
quark NFPDs by q and Ag respectively. Since the cross section is quadratic 
in the parton distributions we find pure gluonic contributions of the form a gg , 
a A g Agi and <Jg& g , similar expressions for the pure quarkonic part, and inter- 
ference contributions a qg , u^qAg, and <J q A g +gAq- Only the interference terms of 
polarized and unpolarized parton densities c g Ag, &gAg, and <J g A q + q Ag depend 
on the alignment of the spins of the incoming electron and proton and thus 
contribute to the spin asymmetries Act := \{<J^i — cr^). (a-^ is the cross section 
for antiparallel lepton spin and nucleon spin.) Comparing the contributions to 
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the spin average we see that a A g A g , CAqAq, and (JAqAg are roughly one order 
of magnitude smaller than a gg , a qq and a qg respectively. 

For our NFPD Models, the gluon contribution dominates the cross section, 
both in the unpolarized and polarized channel. Note that the normalization 
of the cross section depends crucially on the form factors (17). For the used 
input distributions we finally get a spin asymmetry of order 10%. 

The unpolarized and polarized cross sections are plotted in the upper pictures 
in Fig. 6 as functions of the skewedness (. The cross sections receive their 
dominating contribution from the small values of (. In our analysis we restrict 
ourselves on skewedness values larger than ( > 0.1. To go to smaller values of 
£ ultimately would lead to the requirement that Q 2 — > and thus would lead 
to the region of photoproduction. In that kinematical domain the scattered 
electron is not detectable and thus an unambiguous reconstruction of the final 
hadronic state is not possible unless the outgoing proton could be detected 
by a recoil spectrometer. One method to distinguish the elastic events from 
inelastic ones in this case is provided by the existence of a rapidity gap. The 
plots in the middle of Fig. 6 show that the difference between the minimum 
of the meson rapidities and the rapidity of the outgoing proton Ay is for all 
diffractive events larger than 2. In the lower diagrams in Fig. 6 we plotted 
the cross section versus the invariant mass squared of the produced particles. 
These plots show the increase of the ratio between the kaon and pion rates 
with the invariant mass M(h + + h~) as suggested by Fig. 5. The shape of the 
polarized and unpolarized cross sections is similar. 

We have analyzed diffractive 7r + 7r~ and K + K~ production as a tool for the 
determination of polarized and unpolarized nonforward parton distributions. 
We implemented the realistic cuts for the HERMES experiment and required 
sufficiently large pt to justify our leading order treatment. With these cuts the 
cross sections are too small to make this channel a promising one for HERMES. 
Successfull measurements would probably require an ENC-type machine [24]. 
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Fig. 6. The unpolarized (on the left) and polarized (on the right) differential cross 
section for diffractive ir + ir~ (solid circles) and K + K~ (open circles) production 
with respect to the variables C, Ay := min{7/ par ti c i c i,2/ pa rticie 2} - V outg. proton, and 
M 2 {h+ + h~). 
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